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Abstract 

Using the method introduced by Grisaru et al., boundary S matrices for the physical 
excitations of the open Hubbard chain with boundary fields are studied. In contrast 
to the open supersymmetric t-J model, the boundary S matrix for the charge exci- 
tations depend on the boundary fields though the boundary fields do not break the 
spin-S l f/(2) symmetry. 
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Recently, one-dimensional integrable models with boundaries have attracted re- 
newed interest. Those models provide relevant informations for the boundary effects 
on the one- dimensional strongly correlated systems. Among others, as for the bulk 
case, the one-dimensional Hubbard model with open boundary conditions (open 
Hubbard chain) plays an important role in this field. 

In this letter, using the method introduced by Grisaru et al. l| (see also ref. 
0), we study the boundary S matrix for quasiparticles of the open Hubbard chain 
with boundary fields. 

Let us first recall the known facts about the (open) Hubbard chain. The Hamil- 

r(oper 
'bulk 



tonian of the open Hubbard chain with boundary fields is given by = + 



^boundary, where 

#bX en) = "EE (iU + ^Wi.) + U £K- T - l/2)(n n - 1/2),(1) 

^boundary = "M^lt ± n ll) ~ ^h{n L \ ± n L} ). (2) 

Here U is the coupling constant, hi is the boundary field at site I G {1, L}, ipj a (resp. 
ipj a ) denotes the annihilation (resp. the creation) operator of an electron with spin 
° £ {t> 1} a t site j G {1, 2, • • • , L}, and rij a = tpj^ja is the number operator. 

It is well known that the bulk Hamiltonian H^£ a ' (on the bipartite lattice, i.e., 
with even L) possesses an SO(4) = (SU(2) x SU(2))/Z 2 symmetry [[J (see also ref. 
0]). That is, together with the ordinary spin-SU(2) symmetry which corresponds to 
the spin degrees of freedom, is also invariant under the action of the so-called 

7]-SU(2) algebra which pertains to the charge degrees of freedom. The boundary 
Hamiltonian ^boundary ( res P- -^boundary) breaks the r]-SU(2) symmetry (resp. the 
spm-SU(2) symmetry) down to U(l). Quasiperticle spectra of the attractive Hub- 
bard model and those of the repulsive Hubbard model are related by an interchange 
of the spin and charge degrees of freedom ||. Then, in what follows, we restrict 
attention to the Hamiltonian with U > (repulsive case). 

The open Hubbard chain with boundary fields has been solved by the (coor- 
dinate) Bethe ansatz method || |7], ||. The Bethe ansatz for this model provides 
eigenstates of the Hamiltonian which are parameterized by the two sets of 

roots ('rapidities') {kj}^ =1 and {A 7 }^1 1 . Here N is the number of electrons and M 
is the number of electrons with down spin. These roots are subject to the (nested) 
Bethe ansatz equations, 

^(L+Dof, h so (h h x _ n ( A * - sin k J ~ ^/ 2 )( A * + sin kj + ic/2) 

(3(k j: hmkj, h L ) - [[ {As _ gin - - . c/2) {As - sin kj _ . c/2) , (3) 



2 



ji (A 7 -As- zc)(A 7 +A s -ic) (A 7 - sin % - zc/2)(A 7 + sin fc, - ic/2) 

8=1 ( A 7 - A <5 + ic ) ( A 7 + A <5 + ic ) j=i ( A 7 - sin kj + ic/2) (A 7 + sin kj + ic/2) ' 

where j — 1, • • • , N, 7 = 1, • • • , M, and 

c = Z7/2, (5) 
/3(x,h) = - r— . (6) 

Note that, in this model, the solutions of the Bethe ansatz equations are restricted 
as Re(kj), Re(A 7 ) > and kj, A 7 7^ 0. The energy of the model is represented as 

N 

E N = -2 cos kj. (7) 

Next, we shall briefly review the work of Grisaru et al. 0. In 0, Korepin gave 
a general method for exactly extracting the bulk S matrix from the Bethe ansatz 
equations. Then, generalizing this method, Grisaru et al. proposed the method for 
determining the boundary S matrix from the Bethe ansatz equations, and applied 
this method to the open Heisenberg chain with boundary magnetic fields [|T|. Also, 
using this method, Essler et al. calculated the boundary S matrices for the open 
super symmetric t-J model with boundary magnetic fields and those for the open 
super symmetric t-J model with an impurity M. 

An essential ingredient of their method is the following quantization condition 



1 10] for a system of two particles, which has the internal degrees of freedom, with 



factorized scattering on a line of length L; 

e 2ip(ei)L Si2 ( di _ e 2 )K x {e u hi)s 12 (9i + e 2 )K 1 (e u h L ) = 1, (8) 

where 9j is the rapidity of particle j = 1, 2, and p(8) is defined by the expression for 
the momentum of a particle on the corresponding periodic system. Here 5i2(#i — O2) 
is the (bulk) S matrix for the scattering of particles 1 and 2, and K\{Q\,K) is the 
boundary S matrix of the scattering for particle 1 off a boundary with boundary 
field h. Under appropriate conditions on Su(0i — 9 2 ) and Ki(9x,h), the equation 
(H) is equivalent to the following scalar equation (after taking logarithm); 

2Lp(9i) + (bulk two-body phase shifts) 

+ (boundary phase shifts for hi and hz) = (mod 2n). (9) 

Note that, due to the factor Si 2 (9i + 9 2 ) in eq. @, the bulk part of phase shifts 
contains the phase shifts for the scattering of the particle 1 and the mirror image of 
particle 2. 
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On the other hand, if the system is Bethe ansatz solvable, it is possible to derive 
the another condition on p{B\) from the counting function that is defined by the 
Bethe ansatz equations. Then, comparing these two conditions, the boundary phase 
shifts can be evaluated (up to rapidity independent constant) [j]. 

We now turn to consider the boundary scatterings of the open Hubbard chain. 
Since for the open Hubbard chain, it is reasonable to consider the length of the 
system to be L + 1, then we put L = L + 1 in the discussions of the scatterings. 

In this letter, we only consider the case with the bipartite lattice and the half 
filled band, i.e., L even and N = L. In this case, the elementary excitations of 
the periodic Hubbard Hamiltonian transform in the fundamental representations of 
5*0(4) inj. These elementary excitations are called spinons which carry spin 
but no charge and holons/antiholons which carry charge but no spin [|, |i~lf . 
The excitation spectrum can be determined by the scattering of these elementary 
excitations. In ref. || Hj , the bulk S matrix for the periodic Hubbard chain has 



been determined by using Korepin's method. This S matrix has the block diagonal 
form with respect to the scattering of the spin excitations on the spin excitations, 
the spin excitations on the charge excitations, the charge excitations on the spin 
excitations, and the charge excitations on the charge excitations. 

For the open Hubbard chain, the bulk part of the Hamiltonian is also SO (4) 
invariant. Thus, the elementary excitations are still spinons and holons/antiholons. 
However, in our choice of the Hamiltonian, the r)-SU(2) symmetry is broken down 
to U(l). Thus the total 77-spin is not a good quantum number. The boundary S ma- 
trices i^ spin (A, h) and K chSbTge (k, h) for spin and charge excitations, respectively, have 
the following diagonal form, since the Hamiltonian has U(l) x U(l) symmetry 
which corresponds to the preservation of spinon and holon/antiholon numbers; 

K^(A,h) = B(A ° ft) ), (10) 

iWM) = ( c( * ' !) v( l h) ). (ID 

Since the boundary Hamiltonian (|2|) does not break the spin-S'?7(2) symmetry, we 
expect that the boundary S matrix for the spin excitations is proportional to the 
identity matrix, i.e., A(A) = £>(A). In fact, we will confirm this fact. Also we 
define the corresponding boundary phase shifts by the formulae; -4(A, h) = e m( - A ' h \ 
B(A, h) = e tb( ~ A ' h \ C(k, h) = e lc( ~ k ' h \ and V(k, h) = e ld ^ k,h \ From the same argument 
as was given by Grisaru et al. JT|, to determine the above four components, it is 
sufficient to analyze the highest weight states and the lowest weight states of the 
spin (resp. charge) excitation with S = 1 (resp. 77 = 1). Here S (resp. rf) denotes 
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the total spin (resp. 77-spin) quantum number. Note that hereafter we call the states 
which become 77 = 1 states when the boundary fields vanish, 77 = 1 states. Notice 
also that, to study the scattering, we can restrict attention to the states near the 
ground state, i.e., the states which have the microscopic number of holes in the real 
roots. 

Let us introduce counting functions for roots {%} and {A 7 } || [7], ^J. As men- 
tioned above, for later purpose, we only need the real solutions of the Bethe ansatz 
equations @ and (§]). In this case, taking the logarithm of eq. (||) and (H), we have 

Uj = z c (kj), (12) 
A = z s (AJ, (13) 

where z c (k) and z s (A) are counting functions for roots {kj} and {A 7 }, respectively; 
z c (k) = ^-\2kL + -\n(3(k,h 1 ) + -\n(3(k,h L ) 

M 

- 0(2 sin k - 2A S ) + 6(2 sin , (14) 

S=-M 

1 N M 

z s {A) = _[_ £ e(2A-2sin%)+ ]T 6(A - A s ) - 6(A)] , (15) 

Zn j=-N S=-M 



with Q(x) = — 2 tan _1 (x/c). In the above expressions, we have used the 'doubling 
trick', that is, we have put A_,5 = —As,k_j = —kj, and A , k = 0. The two 
sequences of quantum numbers {rij}^ =1 and {J 7 } 7 / £ 1 (we call n-sequence and I- 
sequence respectively) take values in integers, and label the state of the model. 
Remark that n/s, which are defined modulo 2L, take values in < rij < N. Also 
remark that, from the formula |0(x)| < n, J 7 's are restricted as < J 7 < — M(= 
/ max ). For instance, the ground state is characterized by M = N/2 (spin singlet) 
and the configuration rij = j, J 7 = 7. 

We shall also introduce the densities of roots and holes. The number of allowed 
solutions for the Bethe ansatz equations (|3D and ([|) in the intervals (k, k + dk) 
and (A, A + dA) are expressed as L[p{k) + p h {k)]dk and L[a(A) + a h (A)]dA. Here 
p(k) and cr(A) are the densities of roots (filled solutions), and p h (k) and a^A) are 
the densities of holes (unfilled solutions). These are determined by the counting 
functions as follows; 

L[p{k) +p h (k)} = dz c {k)/dk, (16) 
L[a{A)+a h {A)] = dz s {A)/dA. (17) 

Conversely, we can determine the counting functions from the integration of the 
above formulae, if we know the explicit form of p(k), p h (k),a(A), and a h (A). 
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In the thermodynamic limit [L — > oo with N/L and M/L fixed), we obtain the 
following formulae; 

1 r B 

p(k)+p h (k) = - + 2cos£;/ dAa(A)K(2 sin k — 2A) 

7T J-B 

+J-[ T (k, h x ) + r(A;, - ^^K{2 sin fc), (18) 

I/7T L 

<r(A) + a h (A) = 2 [ Q dkp(k)K(2A - 2 sin k) - 1° dA'a(A')K(A - A') 

J-Q J-B 

+lK(A), (19) 

where K(x) = c/[tt(x 2 + c 2 )] and r(x, h) = (h cosx — h 2 ) /(l — 2h cosx + h 2 ). Here 
the charge and spin pseudo Fermi- momenta Q and B, respectively, are determined 
by the conditions 

[ Q dkp(k) = (2N+1)/L, (20) 

J-Q 
rB 

/ dAo-(A) = (2M+1)/L. (21) 

J-B 

Since we have to determine the densities of order 1/L, we may expand p(k) and 
a (A) as 

p(k) = p (k)+ Pl (k)/L + O(l/L 2 ), (22) 
<r(A) = a (A) + a 1 (A)/L + O(l/L 2 ). (23) 

For example, we can easily derive the ground state densities p*o\k) and Oq (A) of 
order O(L ). In the ground state there are no holes. Also, in the half filling case, 
we see Q = tt and B = oo for order 0(L°). Then we can solve eqs. ([18]) and (|19|); 

{g) 1 { cos A; ^j (p) e - J P sinfc - c lpl/2 



/TO = - + rfp — , (24) 

7T 27T J-oo COSh(cp/2) 

( 9)m 1 r ^^o(p)e- ipA 



o-W(A) = — / rfp uv f ; r , (25) 
v ; 2W-00 cosh(cp/2)' v ; 

where Jo(p) is the zeroth-order Bessel function. 

To determine the boundary S matrices for the open Hubbard chain, we now 
proceed to study the excitations which are classified by S z and rj z eigenvalues. Here 
S z (resp. rj z ) denotes the total z-component of spins (resp. //-spins). We will consider 
the spin and charge excitations separately. 
Spin excitation 

We consider the spin excitations. 
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We first investigate the state with S = 1, S z = 1. From this excitation, we can 
determine the component A(A, h) of the boundary S matrix i^ spin (A, h). 

The S = 1, S z = 1 state is obtained by M — > N/2 — 1 with N fixed. In this 
case J max = N — M = N/2 + 1. Thus there are two holes if, 1^ in the /-sequence, 
and the n-sequence do not change. We denote the corresponding spin rapidities A^ 
[a = 1, 2), that is, l£ = z s (A^). The hole densities are thus given by 

p\k) = 0, (26) 

a h (A) = i[<J(A-Aj)+«J(A + Af)+«J(A-A*)+(J(A + A*)]. (27) 
L 

Then we obtain the integral equations for pairs (p (k), cr (A)) and (pi(k), 01(A)) 
with integration boundaries Q and I? which are defined by eqs. (p0|), fl2~T|). Since 
we determine the densities of order 0(1/L), the shifts of the integration boundaries 
from the ground state must be examined of order 0(1/ L). Following refs. Jl0|, [l], ||, 
we assume that, in the thermodynamic limit, the shifts of the integration boundaries 
are of order 0(1/ L n ), (n > 2), as far as the boundary phase shifts are concerned. 
Under this assumption, integral equations can be solved. We then obtain po(k) = 
p ( j\k),a (A) = 4 9 \A),&nd 

, , . cos k [ 1 1 



2c 1 cosh[7r(sin — A^/c] cosh[7r(sin/c + A^)/c] 

1 1 
+ ; f ; . ; ttt^t + 



cosh[7r(sin /c — A^/c] cosh[7r(sin + A^/c] J 

COS^ Z" 00 /" 7r p -2ip(sinfc— sinfc')— c[p| 

+ T^J oo dp j_/ k '\^ k 'M + r(k'M)}- 



COS A; f 00 e -«psinfc-c|p|/2 



cosh (cp) 



2tt 7-oo l + e c W 
+7?-[r(A;,/ii) +T(k,h L )) - 2 cos kK (2 sink), (28) 

27T 

1 /-oo p-«pA 



a x (A) = — / rfp[cos(pA?) + cos(pA^)] 

7T J-00 



1 + e-*! 



1 r djfc _r(A;, /ii) + r(fc, /i L ) 



47rc 7-tt cosh[7r(sin fc — A)/c] 
+ TT- / dp——j-r. (29) 



2tt 7-oo 1 + e- c lPl 

From the above formulae, we can obtain the following equation for the counting 
function z s (A) in the thermodynamic limit; 

-2nz s (A h 1 ) = 2L Ps (A h 1 )+Af 1 (A h 1 )+M 2 (A h 1 ,A h 2 ) = (mod 2vr), (30) 
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where p s (A^) is defined by the expression for the spinon momentum of the corre- 
sponding periodic system |L2|, [|, pTT| ; 



p.(Af) = - r ^ y> sm(M?). (si) 

Jo p cosh (cp 2) 



'o p cosh(cp/2) 
Terms yVi(Ai) and jV 2 {Ai, A£) in (0) are given by 

M(Af) = 7(-2Af/c) + 7(-A^/c) 
1 



dfc[T(fc,fci) + r(fc,fc L )] 

Z7T J-tt 

x<£(-i(Aj - sin fc)/c, 1/4 + i(Aj - sin k)/2c), (32) 
M(Ai, Ag) = 7(-(A?-A^/ C )+ 7 (-(A? + A^)/ C ), (33) 



where 



(f)(x,y) 



00 du(l- e- 2xuJ )e~ 2yul 



uj 



r(i+»)r(y+i/2)' ^ ' 

7(2;) = — 0(ix, (1 — ix)/2) 

r((i-»x)/2)r(i + »x/2) 

~ r((l+zx)/2)r(l-za;/2)- 1 J 

We see that A/^A^A^) are the bulk phase shifts due to the scatterings of the 
particle 1 and 2 , and also the particle 1 and the mirror image of the particle 2. 
Similarly, we can conclude that AZi(Aj') is the sum of boundary phase shifts for the 
scattering of particle 1 off boundaries with boundary fields hi and h 2 . That is, 
j\fi(A) = a(A,hi) + a(A,h 2 ). Therefore we determine A(A,h) up to the rapidity 
independent constant. 

To calculate the component B(A,h) in the equation ([10|), we next consider the 
state with S = 1, S z = —1. We find that, for this spin-SU(2) invariant case, the 
Bethe ansatz equations and energy spectrum of the S — 1, S z — — 1 state are trivially 
same as those for the S — 1, S z — 1 state. Thus we have «4(A, h) = B(A, h). 
Charge excitation 

Next we consider the charge excitations. To determine the two component of the 
boundary S matrix for the charge excitations, we must consider the t] = l,r] z = 1 
state and the 77 = l,r] z = — 1 state by the Bethe ansatz. 

The T] = 1, r] z = — 1 state is obtained by removing two fc's from the ground state, 
i.e., N = L — 2 and M = N/2. In this case, we have 

a h {A) = 0, (36) 
p h (k) = i[S(k - jfef) + 5(k + ifef) + 5(k - k$) + 6(k + *£)]. (37) 
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Similar to the case of the spin excitation, under the assumption for the integra- 
tion boundaries, we obtain the densities 00(A), po(k), 01(A) and p%(k). Results are 
Po(k) = pi g \k), oo(A) = a { 9) (A), and 

~ cos /c f°° g — c IpI/^ 

Pi(k) = ~Lp h (k) — / dp[cos(psin k%) + cos(psin k%)] 

ZTT J-00 



cosh(cp/2) 

rrmk C°° r n p -ip(sink-smk')-c\p\ 

1 ' dp dk'{T(k'M)+T(k'M)\ 



07T J —CO J— 7T cosh(cp/2) 
COS A; f 00 e -«psinfc-c|p|/2 

+^ — / dp- 



2tt 7-oo 1 + e c lpl 
+^-[r(k,h 1 )+r(k,h L )], (3* 

If 1 1 



2c [ cosh[7r(A — sin kf[)/c\ cosh[7r(A + sin ki)/c] 

1 1 

+ . rrrvr + 



cosh[7r(A — sin k^/c] cosh[7r(A + sin&D/e] 
| 1 r dfc _ r ( fc > ^1) + r (^ M 



47rc J-tt cosh [7r (sin fc — A)/c] 

1 poo e -ipA 

+7T-/ dPT-T3CT- ( 39 ) 



2tt i-oo 1 + e c N' 
Also, we have the counting function in the thermodynamic limit 

- 1-nz c {k\) = 2Lpl*=~\k h 1 ) + Mi(k^) + M 2 {k h 1 , k\) = (mod 2tt), (40) 

where p2 z=_1 (^i) * s the quasiparticle momentum of the corresponding periodic sys- 
tern 

pT-\k>{) = -k h x - r^ ^^ ^sin*?), (41) 

jo p cosh (cp/ 2 J 

and 

A^i(^) = 7 (-2sinA; 1 Vc)-2A;f-i[ln/5(A;f,/i 1 )+ln/5(fc 1 h ,/i L )] 
-0(i sin fcf/c, 3/4 - z sin fcj/c) - 6(2 sin fc£) 
-77- r ^ , [r(A;',/i 1 ) + r(A; / ,/i L )]7(-(sinA; / -sinA;f)/c), (42) 
M 2 {k\,k h 2 ) = 7 (-(sin^-sinA; 2 h )/c) + 7(-(sinA;| t + sin4)/c). (43) 

Then, we can determine the component T>(k, h) from A^i(fc^). 

Finally, we have to determine the remaining component C(k,h) in (11). Let us 
study the 77 = l,r) z — 1 state to determine C(k,h). Since the r] — 1,T] Z — 1 state 
is not the regular Bethe ansatz state [§J, we must take the completely filled state 
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\Q) = Ylf=i V4tV4i|0) as the Bethe ansatz vacuum. The Bethe ansatz state with 
2L — N electrons is thus given as 

N 

\®n) = e ii,-w( n b-'.%)ii^ii fl >' ( 44 ) 

o-i,-,o-jve{T,i} i=1 
where n^'s denote the location of electrons on the chain. It is easy to see that the 
eigenvalue of the Hamiltonian H^ +s> for this state is given by E' N = —En, and the 
Bethe ansatz equations are obtained by taking c — > — c in the equations (|3[) and 
Then the problem reduces to find the eigenstates of the attractive Hubbard model 
with the eigenvalues which are given by changing those signs from the correspond- 
ing eigenvalues for the repulsive case. That is, the ground state configuration of 
rapidities for our model is identical to the highest energy configuration of rapidities 
for the attractive case. This is the configuration that all rapidities are real and 
N = L, M = N/2 fI5| . Therefore, the r\ = l,r] z = 1 state is obtained by removing 
two fc's from the ground state configuration. Repeating the calculation similar to 
the case of the rj = 1, i] z = — 1 state, we have 

2vz c (k h x ) = 2Lp^ =1 {k>t) + M'M) + M'M, fc£) = (mod 2vr), (45) 

where p r /, z=l {k 1 l) is the quasiparticle momentum (note that p]! z=1 is different to 

P? =1 (*f) = ti-r- ^rln si ^sin^), (46) 

JO V COSh [CD 2 



and 



M' x {k h x ) = 1 (-2sink h Jc)+2k h l +-[\nP(k^h 1 )+\nf3(k h l ,h L )} 
-0(isin^/c,3/4-isin^/2c) - 9(2 sin k\) 
_J_ f W dk'lrik^h^+rik^hLM-ismk' -sink^/c), (47) 
M' 2 (k^,k^) = 7(-(sin^-sin^)/c) + 7(-(sinA;| t + sin4)/c). (48) 

As for the case of the r] — l,r) z = — 1 state, we obtain the component C(k,h) 
from A^i(fci). 
Boundary S matrices 

Now let us summarize the results. Up to rapidity-independent phase factors, the 
resulting boundary S matrices are expressed as 

K spin (A,h) = e x ^( I J), (49) 
K charge (k,h) = eWOf ^ Q ^ ^^-i J , (50) 
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where 



2X s (A,h) = 7 (-2A/c)+ 7 (-A/c) 
1 



dkr{k, h)<j>(-i(A - sin k)/c, 1/4 + z(A - sin k)/2c), (51) 

71 J -IT 

2X c (k,h) = 7(-2sinfc/c) 

—<p(i sin k/c, 3/4 — i sin k/2c) — 0(2 sin fc) 

— dk'r(k',hW -(sink' -sink)/ c). (52) 

7T J-7T 

It is noteworthy that, in contrast to the case of open supersymmetric t-J model 
0, the boundary S matrix for the spin excitations depends on the boundary field 
although the boundary field does not break the spin- SU (2) symmetry. 

If the boundary fields vanish, the boundary S matix of the charge excitations 
becomes proportional to the identity matrix as expected. The bulk S matrices for 
the Hubbard chain [|5| fll| and the supersymmetric t-J model [[| have the same form 



as that for the XXX chain. However the boundary S matrix for the open Hubbard 
chain has different form with the one for the open XXX model. Full details and 
applications of our results will be published elsewhere. 

Acknowledgment: I wish to thank Dr. T. Yamamoto for discussions and com- 
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